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Abstract 

We develop a deterministic ab-initio model for the input-output relationship of a multiple- 
input multiple-output (MIMO) wireless channel, starting from the Maxwell equations com- 
bined with Ohm's Law. The main technical tools are scattering and geometric perturbation 
theories. The derived relationship can lead us to a deep understanding of how the propagation 
conditions and the coupling effects between the elements of multiple-element arrays affect the 
properties of a MIMO channel, e.g. its capacity and its number of degrees of freedom. 
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u 1 Introduction 

d 

In "wireless communication, systems with multiple antenna arrays at both ends of the link are 
called MIMO (multiple input multiple output). The main interest is to compute the number of 
bits/s which can be transmitted by these systems. It has been shown that this quantity, called 
channel capacity [H HI [31 IH [S] , is linked with the transfer matrix H whose element Hmn is the 
ratio between the complex current intensity /,„ at the load of the to*'* receiving (RX) antenna 
over the current /„ which feeds the n*'* transmitting (TX) antenna. If M and N are respectively 
the number of RX antennas and TX antennas, then "H is a M x matrix. 
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Even if one can experimentally determine H in a given environment when the antennas are fixed, 
its entries show such large variations when moving the antennas arrays from place to place that 
the knowledge of % for some number of locations seems to be insufhcient for having an idea of the 
capacity everywhere. Fortunately, it seems that the eigenvalue distribution of the matrix %*% 
(here * means adjoint) is more robust with respect to the displacements, but this has to be better 
understood. In order to mimic the variations of the entries of "H, one possibility is to replace % by 
random matrices introducing several distributions for their elements (see [SI [7] ) . In sharp contrast 
with this method, we try to understand the variations through the physical study of the wave 
propagation in a given environment. Thus we employ a deterministic modeling method and our 
aim is to describe % by quantities linked with the physical characteristics of the antennas and 
the environment. But certain stochastic aspects can naturally appear if some parameters of the 
scatterers present in the medium are described by stochastic variables (like position, density etc). 
This work was motivated by the existence of an heuristic formula intensively used by the MIMO 
signal experts. It appeared in the works of B.Fleury et al. see e.g. [51 [HI [TOl E]; see also [T^ . 
This formula describes the transfer matrix in the following way: 



%= ( c{nRfA{nR,nT)c{nT)dnRdnT. (i.i) 



where ^It = {4>TtOt) is a direction of departure from the TX side where (Pt^Ot are respectively 
its azimuth and elevation angles and VLr — {4ir,0r) is a direction of arrival at the RX side, 
c{Q.t) and c{^r) are respectively 2 x N and 2 x M matrices whose elements depend only on the 
corresponding antennas diagrams and the geometry of the arrays, while A(flR, fir) is a 2 x 2 
matrix which depends only on the environment. 

Our aim is to understand this formula from first principles, that is solving the Maxwell equations 
in a typical complex environment formed by the TX and RX antenna arrays surrounded by 
buildings , persons, trees ... The antennas are pieces of metal described by their geometry and by 
their conductivity cr, while the scatterers which are made of dielectric and magnetic materials can 
also be described by their geometry, permittivity, permeability and conductivity. 
To determine Hmn one introduces in the second member of the Maxwell equations a harmonic 
current density Ji„, in a small volume constituted by a part of the wire which connects the 
^th gjgj^g^j source to the antenna. The current intensity crossing a section St of the wire is 
In = Jin(a^) • dSr- If one knows the electric fields E induced by this current density, in 
particular at a section Sr of the wire connecting the m*'' RX antenna to its load , there one can 
calculate the current density using the Ohm law Jout = crE and deduce the current crossing this 
section, I„i — a J^^ Ei(x).dSR. Then Hmn = 7^; note that this formula is true when only the n'th 



emitting antenna is fed with current. For the general case see the discussion around (2.201. 
To tackle these calculations, wave propagation experts try to address the difficulties in the following 
manner. They suppose that they know the electromagnetic fields (EM) fields radiated by the n*'' 
TX antenna in the free space. These fields are considered as incoming fields in a scattering problem 
where the scatterers are constituted by buildings, persons, trees... and they try to calculate the 
total fields produced adding the incoming fields and the scattered fields. This is in fact a difficult 
problem which is only solved analytically in some very simple and unrealistic situations. Practically 
it is solved using the geometrical optics approximation (ray tracing). Finally, the total EM fields 
calculated in this way are considered as incoming fields for the m*'* RX antenna and, once more, 
one is faced with a scattering problem where now the RX antennas play the role of scatterers and 
again one has to determine the new total fields, from which one obtains /,„ . 
One of the purposes of the paper is to justify this procedure. We succeed in proving that when the 
antennas are sufficiently small and spatially separated from the surrounding scatterers, the EM 
fields calculated with the procedure explained above appear as the first term in a series giving 
the exact fields. In Lemmas |3.1| and |3.2| we give an approximation for the total Green function 
which puts into evidence the decoupling between the TX antennas, the RX antennas and the 
surrounding scatterers. 

Our proposed methodology is to use the vector potential formalism to solve the Maxwell equations. 
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Our main mathematical tool is geometric pertm'bation theory for closed but non-self adjoint oper- 
ators. These methods have been developed for apparently unrelated quantum scattering problems 
in [121 [TH [T21 US] and they can be used for describing the propagation of electromagnetic waves 
in frequency domain. 

In Corollary 3.4 we derive a stronger version of the Fleury heuristic formula as a corollary of 
Theorem |3.3[ in the case when the antenna arrays are far away from the surrounding scatterers. 
In particular, we link the heuristic "spread matrix function" A{flii,^T) to the Green function of 
the environment alone. 

In this paper the language is deliberately made free of too technical notions of functional analysis 
and operator theory in order to facilitate the understanding for non-mathematicians. Some notions 
concerning the Limiting Absorbtion Principle for the non self adjoint operators coming from the 
vector potential description of the EM fields are not yet proved, so this paper has to be completed 
in this respect in order to get rigorously the set of frequencies for which the Green operators exist 
when the "energy" tends to a real number, as operators in certain weighted spaces. 
On the other hand, this paper is a starting point to more applied engineering works since now we 
understand the link between the heuristic "spread function" in formula (1.1 1 and the scattering 
transfer operator of the environment. 



1.1 Main technical assumptions 

Throughout this paper the configuration space will be M^. The environment is known, which 
means that we assume the knowledge of the dielectric constant e and the conductivity a. They 
are smooth functions of the position and frequency. There are no magnetic effects, that is ii = fiQ 
is constant. 

Let us assume that we know the current density 3 in (x, t) (i.e. in time domain) at the entrance of 
the transmitting antenna. It is completely determined by the transmission device. 

Assumption 1.1. The current density Jj„ has the following generic properties: 

• it is smooth in the time variable; 

• it is compactly supported in x; 

• Ji„(x, t)e~^"*c?t —: J(x, a;) = J(x, — is a smooth function, compactly supported in 
uj € [—uj2,—UJi] U [wi,aj2]. This equality is a consequence of the reality of 3in(jX-,t). 

There is an inherent physical imprecision on the very notion of "external" source current, since 
this quantity is supposed to be completely determined by the transmitter. It is common to assume 
that the current densities on the wires feeding the transmitting antennas are constant on the disc 
of the corresponding cross sections. But this fact does not affect the conclusions of our paper. 

Assumption 1.2. Further technical assumptions on the environment: 

• e(x, io) — e(x, — cj) and (t(x, uj) = (t(x, —uj) are smooth in x. The case with piecewise constant 
e requires a rather complicated regularization procedure which will be considered elsewhere; 

• the dielectric constant of the air is constant and equals eg; 

6('-)(x) := ^ = 1 + Seri^) + Se^i^) + Sen{^), (1.2) 

where all the relative 5e 's are compactly supported perturbations, with disjoint supports. Here 
T means the transmitting antenna(s), R the receiver(s), and M the scatterer(s). 

• the conductivity can be written 

(t(x) = (Tt(x) + crH(x), (1.3) 
where the a 's are compactly supported on the regions containing the antennas. 
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2 The Maxwell equations in the frequency domain 



In the frequency domain, the electric and magnetic fields E and H satisfy the following system of 
equations: 



V • H = 0; 

V X H = J + (jecj + cr)E; 

V X E = -j^ot^H; 

P = V - (eE); 
-jLop = V ■ (J + o-E). 



(2.1) 
(2.2) 
(2.3) 
(2.4) 
(2.5) 



Here we incorporate Ohm's Law in (2.2 1, in the sense that an electric field E generates a charge- 



current density txE in metals. We stress that we do not work with perfect conducting antennas; 
taking <j to infinity is a singular operation within our formalism, and the study of this limit is a 
very interesting mathematical problem which is left as an open problem. We note that one can 
start with infinite conductivities right from the beginning, but the price is that one needs to deal 
with severely ill-posed inverse problems (see e.g. p71 118j ). 

2.1 The magnetic vector potential method 



Equations (2.4) and (2.5) allow us to eliminate the unknown charge density p: 

V • {J + {jeuj + cr)E} = 0. 



(2.6) 



Equation (2.1) allows us to represent H = V x A, where A is a magnetic vector potential. Using 



this in (2.3) gives V x (E -|- j/iQwA) = 0, thus there exists a scalar potential (p such that: 

E + j/iowA = -V0. (2.7) 



Introducing the magnetic vector potential in (2.2) and using (2.7) we obtain: 



V(V • A) - AA = J + {jew + a){-jpoUjA - V^). 



(2.8) 



Note the important thing that if we have a pair (A, (p) which solves ( 2.8 ), then (2.6 ) is automatically 



satisfied because the left hand side of ( 2.8 ) is divergence free. Thus we have the freedom of choosing 



a (p completely determined by A through a Lorenz-type gauge condition, and then solve (2.8) for 
A. Denoting by A^ such a special solution, we choose: 



1 



jeco + (T 



(V • Al) 



With this choice, (2.8) simplifies to 



-AAl = J - jfioLu{jeuj + (t)Al - (V • Al 



(jeuj + a) 



Introducing the notation 



we can rewrite (2.10) as 



fc^(x) := /zoe(x)w^ - jpoUja{x.), 
-AAi - k^AL + (V • Ai)Vln(fc2) 



Now if we have a solution to (2.12 ), then the fields E and H are given by: 

V- Ai 



E = -JpoujAl 
H = V X Al. 



jeu) + a 



(2.9) 

(2.10) 
(2.11) 
(2.12) 

(2.13) 
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Denoting by fcg fJ-a^o^^ and using e'^'")(x) = e(x)/eo (see (1.2|) we have: 

A:2(x) ^kl + (ef'-^Cx) - 1)^2 - j>o^^t(x) fcg + 5fc'(x), 



(2.14) 



where 6k (x.) = (e'^''^(x) — 1)A;q — j/io^cr(x) is smooth and compactly supported. Note that k 
cannot be zero if e*^''^ ^ 0. 

Define the following first order differential operator W in the following way: 

WA:^ -6k^A + {V ■ A}Vln{k^) (2.15) 

and observe that it is relatively compact with respect to —A on the Hilbert space [L^(E'^)]^. 
Moreover, the operator sum —A + is closed on the Sobolev space [_ff^(IR'^)]'^, and if z S C then 
we can define the inverse (—A + W — z)^^ for at least when |/r7i(z)| is large enough. 



2.2 The transfer matrix 



The solution we are looking for in (2.12) is: 



Al = (-A + W^-fc2_j-o+)-iJi„, 



(2.16) 



where we have to give a rigorous sense to the limiting absorption principle. The problem is far 
from being trivial because the perturbation W is not symmetric; technical details and full proofs 
will be given elsewhere. We only want to mention that under rather general conditions, the limit 



in (2.12) is expected to hold outside a discrete set of frequencies. In order to simplify the notation. 



we stop writing 0+. 

Having determined A/,, we can derive E from (2.13): 



1 



jeuj + a 



V-(-A + iy-fc2)-ij. 



(2.17) 



Let us be more precise and assume that we have > 1 transmitting antennas for which we can 
write: 



N 



T _ 7-(") t(") 

•'m — / , ^in "'in > 



(2.18) 



n=l 



where I^^^ (w) is the current running through the tranverse section of the n-th antenna, and 

is its corresponding normalized current density (i.e. the flux integral of any through the 

transverse section of antenna n equals one) . 

According to (1.3), the receiving region has a conductivity a^. Let us assume that there are 



M > 1 disjoint receiving antennas so that we can write: 

M 

''R 



E(m) 
('I 

m— 1 



(2.19) 



Via Ohm's Law Jout = crflE we determine the current density at the m-th receiving antenna. 
Assume that Sm is the transverse section of the m-th receiving antenna. Then the current induced 
in it will be: 



1^7^ = / 4"^ (x)E(x) -dS^Y. ^-4"^ ' 



a(r)(x){-jMo-(-A + M/-fco')-ij|") 



' v.(-A + w^-fcr^j|;^' 



, ("0 
jew + cr)j 



(2.20) 
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The transfer matrix elements T-Lmn are only frequency dependent and give the current induced 
in the m-th receiving antenna when only the n-th transmitting antenna is fed with current. One 
must note that the above formula takes into consideration all possible couplings. In the rest of 
the paper we will try to reduce the complexity of this formula and to arrive to a simpler and 
practically more convenient expression. 

3 Decoupling the receivers and transmitters from the envi- 
ronment 

It is important to realize that the resolvent (—A + W — fco)~^ contains the whole information 
about both the electric and magnetic fields. We will express this resolvent and thus Hmn in a 
different way, which shows a decoupling between transmitters, receivers and scatterers. 
Let us introduce some notation which would describe a situation in which only the transmitter 
would be present: 

fcf,(x) = /cq + SeT{x)kl - j>owcrT(x) =: + Sk'^{x.), (3.1) 

where Sk^ is again smooth and supported only near the transmitter(s). In a similar way we 

introduce the corresponding quantities for M and R. 

The perturbation corresponding only to the transmitter(s) Wt is: 

WtA := -5k^A + {V ln(4)}{V • A} (3.2) 

and similar objects can be defined for R and M . 
Let us introduce the following operators: 

H = -/:^ + Wt + Wm + Wr, Ht:^-A + Wt, Hm:=-A + Wm, i^i? := -A + I^^. (3.3) 

Here Ht only takes into consideration the perturbation induced by the transmitter(s), Hm does the 
same thing for the environment, and Hj^ for the receiver(s). Assume that both the transmitter (s) 
and receiver(s) are separated from all other scatterers, and from each other (see figure 1 for what 
follows) . 

Figure 1. 
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If ^T{RM) E^re bounded open domains completely containing the supports of WT{RM)^ then we 
quantify this separation by: 

dist[nT, f^j?] = Dtb. > 0, min{dist[f}fl,, f^Af], dist[f]T, ^m]} = Dm > 0, (3.4) 
r := 1 + max{diam(r2T), diam(il7?)}, r < mm{DM /2, Dtr/^}. 

Remember that we are interested in the study of the operator {H — z)^^ for z = kg + jO+. 
Denote by xt, XR ^-nd xm the characteristic functions of fir, ^r and f^M respectively. Define a 
smooth function < Jm < 1 which enters in ^It and il^ and "touches" neither the transmitter (s) 
nor the receiver (s): 

Ja/(x) = 1, x^rfTUl^fl, (3.5) 
Jm Xsupp(Wt) — Jj^t Xsupp{Wr) = 0- (3.6) 

The second condition is possible because the fi's completely contain the two antenna systems. Let 
us also note the identity: 

JaKx){1 - Xt(x) - Xi?(x)} = 1 - Xt(x) - Xi^(x), x G (3.7) 

Since we assumed r > 1, we can choose a smooth function < J < 1 such that 

J(x) = 1 if |x| < r - 1, J(x) if |x| > r. (3.8) 

Take a point in as the origin of coordinates and let e e fi/j. Define 

Jt(x) := J(x), J«(x) := J(x - e). (3.9) 



Our choice of r in (3.4 1 insures that r — 1 is larger than the diameters of both fir and Qji. Thus: 

JtXt = XT, JrXR = XR- (3.10) 

The support of any derivative of Jt is contained in the annulus r — 1 < |x| < while the support 
of any derivative of Jji is contained in the annulus r — 1 < |x — e| < r. 

Let us denote by XTa the characteristic function of the spherical annulus r — 1 < |x| < r, and by 
Xr^ the characteristic function of the spherical annulus r — 1 < |x — e| < r. Then we clearly have 
the identities: 

XtA^)D''Jt{^) = i^"JT(x), xrJ^)D^M^) - ^"Jfl,(x). (3.11) 

The supports of XTa and XRa are disjoint from each other, and are situated in the free space. 
If z has a sufficiently large imaginary part, we can define the following bounded operator (see 
Q): 

Siz) := JriHr - z)-\t + Jm{Hm - ^)"'(1 - Xt - Xr) + Jr{Hr - z)-\r. (3.12) 
We then have: 

{H - z)S{z) = 1 + Kt{z) + Km{z) + Kr{z), (3.13) 

where 

Kt{z) = [-A, Jt]{Ht - z)-\t, Km{z) = [-A, Jm]{Hm - z)-\l - xt - Xr), 

Kn{z) = [-A, Jn]{HR - z)-\r, (3.14) 

and this is because Wtirm) commute with Jt(rm)- For example, 

[Wm, Jm]A = {Vln(4^)}{VJA/ • A} = 
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due to disjoint support properties of Vfc^j and VJm- 
If 13(2;) I is large enough, then one can prove that 



xnso,{\\KTMK, 



Ml 



\Kt. 



< 1/10, 



thus 



{H - = S{z)[l + Kt{z) + Km{z) + Kn{z)]-^ 

= S{z) -{H- z)-'[Kt{z) + Km{z) + Kr{z)]. 



(3.15) 



This equahty can be extended to all z where both sides make sense, that is outside of a discrete 
set of singularities. 

Let us observe an important symmetry property coming from time reversal invariance. Define 



VFA := -V ({Vln(fc2)} . A) -5k^A 



(3.16) 



This operator is depending on frequency through fc, and due to Assumption 1.2 and (2.11 ) we may 
write an important identity for the adjoint: 



If H^^ := — A + W-^, then another consequence is: 

{(i?_^-z)-i}* = {H^-z)-\ 
which is also true for T, R and M alone. 



Write (3.151 with z and — w, then take the adjoint and use (3.181. We obtain: 

{H^ - z)-^ - {S-u{z)Y - {Kt{z, -oj) + Km{z, -uj) + Kr{z, -tu)}*{H^ - z)-\ 
We introduce the notations: 

S{z) Xt(Ht - zy^Jr + (1 - XT - Xr){Hm - z^^Jm + Xr{Hr - z^^Jr, 

and 



(3.17) 
(3.18) 

(3.19) 
(3.20) 



Kriz) 
Km{z) 
Kb{z) 



= {Kt{z,~lo)Y ^Xt{Ht~z)-\^,JtI 

= {Km{z, -uj)}* = (1 - XT - Xr){Hm - z)-^[\ J. 

^ {KR{z,-Lo)y = xr{Hr - z)-i[A, JflJ. 



Then (3.191 can be written in a more compact way: 



(3.21) 
(3.22) 



[H - z)-i = S{z) - [Kt{z) + Km{z) + Kr{z)]{H - z)-\ 
Here is the first of our technical results: 

Lemma 3.1. The following "almost decoupled" formula holds: 

Xsmpx>(Wr){H - z)~'^XT = Xsupx>(Wn)^R{'A ' XrAH - z^^xt^ ■ Kt{z)xt- (3.23) 
Proof. There are several important things to note here. First, we have the support condition: 

supp(J„0 C supp(xt)- (3.24) 

Second, due to the support properties of our various cut-off functions we have — K\,j — Kj^ — 
KtKr = KrKt — 0. Now if we introduce (3.15) in the left hand side of (3.23), we see that the 
term Xs\i-pp{Wb)S(z)xsup'p(Wt) — ^ because of various support properties. We obtain: 

Xsupp(iv„)(^f - zy^XT = ~Xsux>v(Wr){H ~ zy^KT{z)xT- (3.25) 
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Now use (3.22) in the right hand side of (3.25). We obtain: 

Xsupp{Wr){H - zy^XT = -~Xsnpp{WR)S{z)KT{z)xT 

+ Xsnpp{Wn)KR{z){H - z)'^Kt{z)xt. (3.26) 

But the first term on the rhs of the above equahty is zero, again due to the supports. Finally, use 
(3.11) in the above equation, and the proof is over. □ 



Remark. Note that Xsnpp(w b.){H — z) "^xt is the operator entering in (2.20) giving the general 
transfer matrix element, because xtJih = Jin and Xsupp(Wr)'^R = ^R- Although (3.231 is an 



exact formula, it is not very useful yet because in the middle of the right hand side we still 
have the full resolvent and not just the resolvent corresponding to the environment. But in the 
next subsection we will show that XR^iH — z)^^XTa can be better and better approximated with 
Xr^{Hm — z)~^XTa if the space occupied by the antennas become smaller and smaller compared 
to the distances between different objects. 



3.1 The study of 



XT, 



Remember that xt^ and XRa are the characteristic functions of two spherical annuli which are at 
a distance r > 1 from both the transmitter(s) and the receiver(s). If the physical space occupied 
by our antennas becomes very small (mathematically this means that the volume of the supports 
of xt and xr i-6- of ^It and fij^ is very small), then it would be natural to be able to approximate 
Xr^{H — z)~^XTa with XRai^M ~ z)^^XTa- Let us show this here. 

Lemma 3.2. Let Vj. :— Yo\{Qj^) and Vt '■= Vo^ri^). Then outside a discrete set of frequencies, 
the vector potential can be approximated in any C'^(riij) norm with k > in the following way: 

XrA - KR{kl+jO+)xRAHM ~ kl - jO+]-\T^KTikl + jO+)3in - o(max{i;T, «r}). (3.27) 



Proof. In the Appendix we have formulated a Lippmann-Schwinger type representation of the total 
resolvent in the presence of N perturbations, see (5.11). We want to particularize that formula 
for A^ = 2 objects. We put Hq = Hm = -A + Wm, we take Wt to be Wi, and Wr will be W2. 



The formula (5.11) reads as 



R{z) = iHM-zr'~Y.J2(HM 



z) ^XnK 



{z)Xm{HM — z)' 



This implies: 



XRa, [H - z) ^XT, - Xfl, [Hm - z) ^xt. 



(3.28) 



(3.29) 



T. 

m—1 n— 1 



i{z)Xni{HM ~ z) ^XT„ 



Now the idea is to show that the right hand side of ( 3.29 ) is small when the volume of the supports 



of xi and X2 are smaller and smaller. We need to estimate the norm of the operator 

XRai^M - zy^XnAnm{z)Xm{HM - zY^XT^- 

The factor XrA^m — z)^^Xn is a Hilbert-Schmidt operator and its Hilbert-Schmidt norm tends 



to zero like 



A similar result holds true for Xm{HM — z) Xt„- It means that the right hand 



side of (3.29) is close to zero when the linear dimensions of our antennas become very small. 



We can thus write: 

XRAH-{kl + ]Q+)y 



Vt„ = XrAHm - {kl+3<d+)) + O (max{Vt^I, v^}) 



(3.30) 



outside of a discrete set of frequencies. This proves the lemma for the norm; one can actually 
show that this convergence is also true in any C*^ norm; the ingredients are the elliptic regularity 
and the fact that the supports of xt, X-Rj XRa and XTa are disjoint. The proof is over. □ 
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Remark. An obvious interpretation of this formula is the following: the input current is 
transformed into a signal by Kj- (only depending on the transmitter(s)) and sent into the annulus 
given by xt^,- From there, Hm scatters the signal into the observation region of the receiver(s), 
or XRa- Finally, Kr takes over the signal and sends it to the receiver (s). Note that the diameters 
of xt„ and XRa do not have to be large, and this is exactly what happens when some scatterers 
are close to our antennas. But the linear dimensions of the emitting and receveing antennas 
have to be small in order to be sure that we can approximate XRai^ ^ i^o +iO+))^^XTa with 
XrAHm - (kl + jO+))-^XT^- 

Another important observation: if the linear dimensions of the antennas are important relatively 
to the other distances in our decomposition, then the expression in (3.27) is not correct. We would 
need to take into consideration the complete formula for the full resolvent (3.29), because we can 
no longer ignore the higher order coupling between the emitting and receiving antennas given by 
the exact formula (3.29). 



3.2 The main theorem 



Remember that the transfer matrix in (2.20) is completely characterized by the value of A = 
(ff - fcg - jO+)-iJ 



In this subsection we assume that the antennas are small enough so that 
(3.27) makes a good approximation. In order to simplify notation, we write zq instead of + jO+. 
All the resolvents we have considered until now are integral operators in the following sense: 

3 

{(iJ-^o)"'«'}a(x) = V / G,t(x,x';zo)«'t(x')dx', se{l,2,3}. (3.31) 

Their integral kernels are smooth functions of x and x' outside the diagonal x = x', due to general 
elliptic regularity results. We will now express the quantity in (3.27) with the help of the various 
integral kernels appearing in that equality. Using (3.21) and (3.14) we may write (here x e flu) : 

dJ. 



E 

t2 



d dJR 

dus^ dusi 



d 



dus^ dus^ 



E 

t3 



E 

S2 



d dJr 

dvso dvs^ 



OJt d 



V / dxTG'lft^Jv,y;zo)Jin,t4(y)- 



(3.32) 



In order to simplify the notation, we will assume that the transmitters are located near the origin 
of coordinates, while the receivers are located near e. In this formula we will choose Jt to be a 
moUifier of the characteristic function of the ball Br{0) = {|x| < r}, and Jji to be a moUifier of 
the characteristic function of Br{e) — {|x — e| < r}. In Br{0) we may choose to work with local 
spherical coordinates {p,xt), with p > and x^ G 5*^. The same thing can be done for Br{e), 
and denote its local spherical coordinates by {p',Scfj). We choose Jt and Jji to be radial in these 
coordinates. Then as they converge towards the characteristic functions, one can prove that the 
expression in (3.32) converges to: 



^^E 

t2 



din I dp'Gifl (x, e + p'in; ^o) | ^<5(p' - r) + Sip' -r)^ 



V / dx / dpG\ftl{e + p'±R,p±T\zo) 
^Js^ Jo 



^^S{p-r) + Sip-r)^^ 



E 



dyGlJApx,y; zo)Jin,tAy)- 



(3.33) 
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Perforin the radial integrals and write At-^ (x) as a sum of four terms: 
At^ (x) « -r"^ y2 / d^rdxR 



- {dp'G^fl e + p'ifl; zo)}p'=rG|f4^ (e + ri^; rir; zo){i9pG^ft^, (pxr; y; zo)}p=r 
+ G|ft|, (x, e + rifl,; 2:o){9p'G|ff3' (e + p'±R, rxr; 2;o)}p'=r{9pG|ft^^ {pxr; y; zo)}p=r 
-G^ft2(x; e + rxn; zo){dl, pG[^2{e + p'Slr, p±T] zo)}p=p'=rG^^t'jrxT; y; ^o)} 
Jin,u(y)- (3.34) 
Let us introduce the notation (boldface G indicates 3x3 matrices): 

c ■^)-f {5pG(^)(pxT,y;zo)}p=. \ ..._f {9p'G(«)(x, e + p'i^^; zo)}p'=. A 

V G^'^(rxT;y;zo) / V ' i^; e + txr; zq) J 

(3.35) 

and 

r G(Af)(e + rXfl,rxT;zo) -{apG(*^)(e + rx^, pxy; zo)}p=. A 

yvnx^,XTj. _{a^,GW(e + p'xfl,rxT;zo)}p'=. {9^,.pGW(e + p'X^; pir; ^o) J " 

(3.36) 



Then the key equation (3.271 can be rewritten as 



A(x) = {(ff-zo)-iJi„}(x) 

sa r* / dy / / dxTrfx7j{cfl(x77; x)}*A^(xfl,; Xt)ct(xt; y)Jin(y), (3.37) 

where the transposition operation is considered with respect to the 2x2 structure. 

Already here we can see the complete separation between the transmitters, receivers, and the 



rest of the scatterers in the environment. Introducing this formula back into (2.201, we prove the 
following theorem: 

Theorem 3.3. In the case in which the linear dimensions of the emitting and receiving antennas 
are small, we can approximate the transfer matrix elements by the formula: 

■Hmn^r^ ! [ dXTdx^(g^"'(xi?),A^(i^;xT)g^"^(xT)), (3.38) 
where (•, •) denotes the usual dot-product in M.^ , while 

g(?)(xT) := J CT(xT;y)j|:^(y)dy, (3.39) 
g(,"'(iH) := ^ dSra ■ |-jpoc^{ch(x«;x)}* + V l-^V • {cfl(i;^; x)}*^ | ^"^(x) 

are two six- dimensional vectors only depending on the transmitters and receivers respectively. They 
characterize the radiation pattern of our transmitting and receiving antennas, and they do not 
change with the environment. 

3.3 Application: the spread function in the case of distant scatterers 

The key mathematical object which completely characterizes how the various scatterers affect 
the signal is the Green function of the environment G'^^^(x;x'; fcg). We will try to obtain some 
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relatively simple formulas for this Green function. We have already given in the Appendix a fairly 
general expression for the resolvent in the presence of N scatterers, see (5.10) and (5.111. 

We would like to give an even simpler, yet generic expression of this Green function when the 
distance between the scatterers and the two balls Br{0) and Br{e) is large enough. As before, 
let us use zq = fc§. Assume that the total number of these scatterers is > 1. The free Green 
function is given by the well known formula 



G^;-(x;y;fco^) = 5, 



pifeo|x-y| 



^''47r|x- 

Let us denote by Ro{zo) its associated operator. 
Then (5.11 1 gives: 

s s 

{Hm - zq)~^ 



H 
c 



(3.40) 



^o(^o) - 2^2^Ro{zo)XuAuv{zq)XvRq{zo) 

U—1 V—1 



Rq{zo) - Ro{zo)TMizo)Ro{zo), 



(3.41) 



Here the operator Tm — XuAuvXv contains all possible interactions between various 

scatterers present in the environment. 

Now assume that the distance between scatterers and the observation points rx^ and e + rxjj 
is much larger than the radius r (see figure 2). 



Figure 2. 




Then if y and y' are in the support of the scatterers we may write: 

|rxT - y'l = |y'| - ricT • y' + OdyT'), y' = ^ 

and 

«' + rxB. - y| = |y - e| - rxji. • (y - e) + 0(|y - e|"^). 

12 



Therefore 



Gf-°(e + Tin; y; kl) = <5,, 4n\y-e\ + ^^'^ " ^^l"'^' 

pj*:o|y'L-j'£o'-XT-y' 

Gf™°(y'; ricT; A:^) = ,5.,^ + 0{\y'\-'). (3.42) 

Using some elliptic regularity estimates and integration by parts, one can prove the existence 
of some 3x3 matrices a-uviymyv', zq) and b„(y„; zq) which are jointly continuous as functions of 
y„ and y^, on the support of Xu and Xv such that: 

{Ro{zo)XuAuv{zQ)xvRQ{zQ)}{e + rx_R,rxT) 

r r g-j*:o''Xn-(y„-e)-jfcorXT-yi. 

za) + b(y„; ZQ)6{yu - yv)} dyudy^ 

Jsupp(x„) ^supp(x„) I67r^|y„ - e| |y.„| 
+ 0{l/Dl,), (3.43) 



where Dm is as in (3.4 ) (i.e. the minimal distance between the emittors/receivers and the scatterers 
of the environment). The matrices a.uv{-, ■', zq) and hu{-', zq) contain the full scattering information 
and depend on all scatterers not just on u and v. Note though that the dependence of x and x' 
is now very explicit. 



Using (3.43) in (3.41) allows us to introduce a single 3x3 matrix kernel tA/(y,y';zo) which 



contains the full scattering information of the medium: 

G(*^)(e + rxfl,rxT; k^) = Gf-°(e + rx^^rxT; k^) (3.44) 
-/ / T^r^ ryvi tM{y,y';zo)dydy' + 0{1/Dl,), 

where xm is the characteristic function of all scatterers, while tM{y,y'',ZQ) has the structure: 

tj\/(y, y'; zo) = ^(y, y'; Zq) + 0{y; ZQ)6{y - y') 

where t(-, •; zq) and 6{-; zq) are continuous on their domain of definition. These functions play the 
role of a scattering kernel, and contain the full scattering information of the medium. 



When we introduce this formula in (3.36) we obtain two contributions: one from the "empty" 
space and another one coming from the scatterers. The same happens with the transfer matrix. 
This scattering generated contribution can be expressed as: 

K^T^r^l f d^Td^n f f dydy'(g^^\±R),M{y,y';±R;±T)g^^\^T)) 

J J Jsupp{xm) Jsupp{xm) 

(3.45) 



where 7M(y, y'; x^j; x^) is constructed by introducing the scattering contribution from ( |3.44 ) in 



(|3.36|) which gives: 

jrfeoXn-(y^)-irfeoXT'y' 



e 



M{y,y';iR:iT) ■■= ,1 11 n 

167r2|y- e| |y'| 

tM(y^y';zo) jfco(xT ■ y')tM(>sy'; ^o) \ ^^^^^ 

jko±R ■ (y^)tM(y,y'; zq) -kK-kx ■ y')(fcoXfl • (y^))tM(y,y'; zq) ' ' 

We are now ready to formulate the main result of this section: 
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Corollary 3.4. The scattering contribution H^nn^ transfer matrix element can be expressed 

as 

nTn'= [ [ dnRdnT(h'^R\nR),A{nR,nTW^\nT)) + oii/Dl,), (3.47) 

where (Q,t) *s « six dimensional vector which can be interpreted as the signal sent by the 
transmiter in the direction VIt, then A{Qr, fir) is a 6 x 6 spread matrix only depending on the 
scatterers, and finally h^^(J7^) is a six dimensional vector describing what the receiver got from 
the direction D,r. 



Proof. Going back to (3.45) we may assume that: 

g^i^r) = ( p^fgj ) , S^rH^r) = ( ilii^n) H2(iK) ) (3.48) 
where the F's and H's are some three dimensional vectors. Then we have: 



Sr 'iiR),M{yy;±R;±T)g>r'i^T)/ - i67r2|y-e| |y'| 

X {Hi(xfl) •tM(y,y';2;o)Fi(xT) + Hi(xi^) • tM(y,y';2o)F2(xT)jfco(xT ■ y') 

+ jko{xR ■ (y^))H2(xfl,) • tM(y,y';2o)Fi(xT) 

- fc2(x3, . p)(ifl • (Sr^))H2(xK) • tM(y,y'; zo)F2(xt)}. (3.49) 



Now we integrate with respect to the angles xr and xt in (3.45). It is useful to define local 
spherical coordinates near both the transmitters and the receivers. Introduce Qt y' and 
st ■= \y'\, and Hr :— y — e and sr := |y — e|. Now we can define: 



( Is^ e"^'^^o*-^-Fi(xT)dxT , , . 



Am) 



(fit 



Jg, e-J"'-'^"*«-"«)Hi(xfl)dxfl, /g, e~^-'^oScn-nn)jko{xR ■ f]fl)H2(ifl)difl, 



(3.51) 



and 

AinR,nT) 



dsj 



ds' 



r'^SRST 
167r2 



^m{sr, flR] ST, ^t; zo) tM{sR, ^r; St, SIt; zo) 
^MisR, ^r; St, f^r; ^o) tMisR, ^Ir; st, ^t; zq) 



(3.52) 



Then (3.47) follows immediately. Note that the integrand in the double integral defining the 
spread matrix A{^r, ^It) is different from zero only on the compact radial supports of the total 
scatterer defined by the intersection of the scatterer with the direction ^It seen by an observer 
placed at the origin and by the intersection of the scatterer with the direction Qr seen by an 
observer placed at e. 

□ 



4 Conclusions and open problems 

1. Starting from the Maxwell equations and the Ohm's Law we give a rigorous analysis of the 
input-output relationship of a MIMO system as a direct, well-posed problem. The main 
observation is that we can decouple the group of transmitters from the group of receivers 



and scatterers if the linear dimensions of all our antennas are small enough, see (3.27) and 



14 



the discussion around it. Even if this smallness condition would not be satisfied, we could 
in theory give higher order corrections with respect to the decoupled case. 

We stress that the transmitting (receiving) antennas are allowed to interfere among them- 
selves, so in principle the coupling in between the transmitting (receiving) antennas is taken 
into account. 

In the decoupled case, we can analyze the transfer matrix and identify in it the spread 
kernel due to the environment alone. The most important result of our paper is contained 



in the equations (3.38) and (3.391; there we do not need the scatterers to be located in the 
far-field region of the receiving and transmitting antennas. In the particular case in which 
the scatterers are far away, then our formulas simplify and we can recover previously known 



empirical results derived in [51 [ini El [3 [H] ; in this case, our results are given in (3.50)-(3.52). 

In a future work we will investigate the behavior of the spread kernel as a function of 
angles and frequency for scatterers which are not necessarily far away from the emitting and 
receiving antennas. 

3. Our formalism does not (yet) allow ideal metals {a = oo) or discontinuities in e's. Investi- 
gating how our formalism behaves when one takes the limit of non-smooth coefficients is a 
very interesting and difficult problem of operator theory and functional analysis, which will 
be investigated elsewhere. 

4. We have not elaborated on the question of the limiting absorption principle formulated in 



(2.16); complete proofs based on the analytic Fredholm alternative will be given elsewhere. 



We think that our new understanding of the construction of the transfer matrix H paves the 
way for the study of the behavior of the capacity in the case the number of antennas grows 
in a definite volume (see [TS] for a discussion on this subject). 



5 Appendix: A Lippmann-Schwinger type equation for the 
resolvent 

Let us consider an operator H = Hq + J2'k=i where Hq is a "nice" elliptic second order 
reference differential operator (the typical example is —A), and the perturbations Wk are first 
order differential operator with smooth coefficients. The supports of the coefiicients of Wk are 
disjoint from those of the coefficients of Wj if fc ^ j. Denote by Xm the characteristic function of 
a ball completely containing the support of the coefficients of Wm- Denote by Dmn the distance 
between the supports of Xm and Xn- 

Denote by R{z) := (H—z)^^ and by Ro{z) := {Hq—z}^^ whenever the two inverses (resolvents) 
exist. 

5.1 The case of just one scatterer 

Assume = 1. Choose z in the resolvent set of Hq. If the imaginary part of z is large enough, 
then ||Wi(i/o — < 1 and {H — z)~^ exists as a bounded operator in (L^(K'^))'^. The second 

resolvent identity reads as: 

R{z) ^ Ro{z) - Rq{z)WiR{z). (5.1) 
Multiply the above equation with Wi at the left, and write 

WiR = WiRo - WiRqWiR. 

Using that Xi^i = W^i we have: 

WiR - (Idi + WiRoXiV^WiRo, (5.2) 
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where Idi is the identity operator in (i^(supp(xi)))'^. The inverse (Idi + W^ii?oXi)~"'^ (if it exists) 
is to be taken in (L^(supp(xi)))'^. We will always assume the existence of this inverse; generically 
this is true if Wi is relatively compact to Rq and one can apply the Fredholm alternative. 

Note that if we know WiR, then we know R everywhere in the space because we can replace 
(5.2 1 in (5.1) and obtain: 

R{z) = Ro{z) - Ro{z)xi{Ui + WMz)xi)-^WMz) 

^ Ro{z) - Ro{z)Ti{z)Ro{z), 
Ti{z) := xi(Idi + WMz)xi)-^Wi. (5.3) 

5.2 The case of several scatterers 



Let N > 2. The equivalent of (5.1 1 reads as: 

JV 



i?(z) = Roiz) - ^ Roiz)WraR{z). (5.4) 



7n—l 

We multiply with Wn at the left on both sides of the above equality and obtain: 

N N 



WnR{z) = WMz) - J2 WMz)W^R{z) = i^"™ - WnRo{z)xm} W^R{z). (5.5) 

m— 1 m— 1 

Denote by Hn ■— ffiI^i(-^^(supp(xA;)))'^. Define the bounded operator M given by 



M.kiz) := W,Ro{z)xk- (5.6) 
Define the immersion operator 



N 



n=l 



Denoting the total perturbation by = X!!^=i the equation (5.5 1 can be seen as: 

JnWR{z) = JnWRo{z) - M{z)JnWR{z). (5.8) 

Denote by Id the identity operator in "Hat. If Id + M{z) is invertible in Hn, then the above 
equation can be rewritten as: 

JnWR{z) = {Id + M{z)y^ JnWRo(z). (5.9) 



Moreover, equation ( |5.4[ ) can be written as: 

R{z) = Ro{z) - Ra{z)J*^JNWR{z) 



^R^{z)~Rfi{z)J*j,{l<l + M{z)} ^JNWRoiz) (5.10) 



where in the second line we used (5.9). It follows that the full resolvent can be always written as: 

N N 

R{z) = Ro{z) Ro{z)XnAnm{z)XrnRo{z), (5.11) 

71—1 m—1 
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where the operators Anm(z) act m (L^(M'^))^. 

Now let us assume that the distance between any two different supports is large. We can 
split the operator in a diagonal and off-diagonal part: 



Md{z) := ®'^=iWnRo{z)Xn, Mo{z) := M{z) - Md{z). 



(5.12) 



Clearly, the norm of A4o{z) becomes smaller and smaller when Dmn becomes larger. For example, 
if Hq — —A we have that ||H/^mJ^o('2^)Xn|| < C/Dmn- Thus if the minimal distance between any 
two scatterers becomes larger than a critical value, we can write: 



{Id + M{z)}-^ - E(-l)' {Id + Md{z)r^ {Mo{z)[ld + Md{zr')y ■ 



j=0 



The operators [Id + Md{z)] ^ are diagonal and given by: 

[Id + Md{z)]-^ = ©^=l[Idn + WnRo{z)Xnr^- 



Making the analogy with (5.3 1 we introduce the notation 

r„(2) := [Id„ + WMz)Xnr^Wn 



(5.13) 



(5.14) 



(5.15) 



where T„ is a bounded operator in (L^(supp(x„)))'^. Then introducing (5.15) and (5.13) in ( 5.10| ) 
we obtain (we drop the z dependence for simplicity): 



N 



~ N 



R — Ro — / , RoTnRo — > ( — 1)/ RoTjnoRoTjn^ RoTrmRo, 

^ — ^ ^ — ^ ^ — 'mo, mi,.. .mi — 1 
n— 1 i—1 



(5.16) 



where the symbol ^ means that the sum is performed on indices which obey mg ^ mi, mi ^ m2, 

mi-i ^ mi. 
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